The magnetic neutron scattering resonance of high-Tc superconductors in external 

magnetic fields: an SO(5) study 
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The magnetic resonance at 41 meV observed in neutron scattering studies of YBa2Cu307 holds 
a key position in the understanding of high-Tc superconductivity. Within the SO (5) model for 
superconductivity and antiferromagnetism, we have calculated the effect of an applied magnetic 
field on the neutron scattering cross-section of the magnetic resonance. In the presence of Abrikosov 
vortices, the neutron scattering cross-section shows clear signatures of not only the fluctuations in 
the superconducting order parameter i/;, but also the modulation of the phase of -i/; due to vortices. 
In reciprocal space we find that i) the scattering amplitude is zero at {n /a,n/a), ii) the resonance 
peak is split into a ring with radius n/d centered at (vr/a, vr/a), d being the vortex lattice constant, 
and consequently, iii) the splitting n/d scales with the magnetic field as y/B. 
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Soon after the discovery of high-Tc superconductivity 
in the doped cuprate compounds, its intimate relation 
to antiferromagnetism was reahzed. A key discovery in 
the unraveling of this relationship was the obsjervation of 
the so called 41 meV magnetic resonanceETH later also 
denoted the tt resonance. In inelastic neutron scatter- 
ing experiments on YBa2Cu30'3-|at temperatures below 
Tc ^ 90 K, Rossat-Mignod et alB found a sharp peak at 
hu! ~ 41meV and q = (7r/a,7r/a), a being the lattice 
constant of the square lattice in the copper-oxide planes. 
Latetits antiferromagnetic origin was confirmed by Mook 
et ala in a polarized neutron scattering experiment and 
subsequently Fong et alB found that the magnetic scat- 
tering appears only in the superconducting state. Re- 
cently, Fong et aln have also observed the n resonance in 
Bi2Sr2CaCu208+5, which means that it is a general fea- 
ture of high-Tc superconductors and not a phenomenon 
restricted to YBa2Cu307. This gives strong experimen- 
tal evidence for the tt resonance being related to an- 
tiferromagnetic fluctuations within the superconducting 
state. Conversely, it may be noted that angular-resolved 
photoemission spectroscopy has shown how the single- 
particle gap within the antiferromagnetic state inheijils 
the d-wave modulation of the superconducting state.ErQ 

A number of differenL models have been proposed to 
explain the tt resonance.0 Ej In particular, Zhang was in- 
spired by the existence of antiferromagnetic fluctuations 
in the superconducting state to suggest a unified SO (5) 
theory of antiferromagnetism and d-wave superconduc- 
tivity in the high-Tc superconductorsJla It is of great in- 
terest to extend the different theoretical explanations to 
make predictions for the behavior of the tt resonance e.g. 
in an applied magnetic field. An experimental test of 
such predictions will put important constraints on theo- 
retical explanations of the tt resonance in particular and 
of high-Tc superconductivity in general. In this paper 
we treat the tt resonance in the presence of an applied 
magnetic field within the SO (5) model. 



Zhang proposed that the cuprates at low temper- 
atures can be understood as a competition between 
d-wave superconductivity and antiferromagnetism of a 
system wliich at higher temperatures possesses SO (5) 
symmetry£3 The SO (5) symmetry group is the min- 
imal group that contains both the gauge group U(l) 
[=S0(2)] which is broken in the superconducting state, 
and the spin rotation group SO (3) which is broken in 
the antiferromagnetic state. Furthermore, the SO (5) 
group also contains rotations of the superspin between 
the antiferromagnetic sector and the superconducting 
sector. The relevant order parameter is a real vector 
n = (711,712,713,714,715) in a five dimensional superspin 
space with a length which is fixed (|n|^ — 1) at low tem- 
peratures. This order parameter is related to the complex 
superconducting order parameter, -0, and the antiferro- 
magnetic order parameter, m, in each copper-oxide plane 
as follows: "0 = /e*"^ = tii -I- in^ and m = (7^2, '^■3, 't-4)- 
Zhang argued how in terms of the five dimensional su- 
perspin space one can construct an effective Lagrangian 
£(n) describing the low energy physics of the t-J limit of 
the Hubbard model. 

Two comments are appropriate here. Firstly, we note 
that relaxing the constraint |n|^ = 1 in the bulk super- 
conducting state will introduce high energy modes, but 
these can safely be ignored at low temperatures. More- 
over, they do not alter the topology of vortices in the 
order parameter, which is our main concern. Secondly, 
one may worry that results obtained from a pure SO (5) 
model deviate substantially from those obtained from 
the recently deyeloped, physically more correct projected 
SO (5) theoryE£l. However, the two models are only signif- 
icantly different close to half filling, and our study con- 
cerns AF-modes in the bulk superconductor in a weak 
magnetic field, a state which although endowed with the 
topology of vortices is far from half filling. For simplic- 
ity, we thus restrict the calculations in this paper to the 
original form of the SO (5) theory. 
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In the superconducting state the SO (5) symmetry is 
spontaneously broken which leads to a "high" energy col- 
lective mode where the approximate S0(5) symmetry al- 
lows for rotations of n between the superconducting and 
the antiferromagnetic phases. These rotations have an 
energy cost Tiw^ corresponding to the tt resonance and 
fluctuations in n will thus give rise to a neutron scattering 
peak at TilOtt which, through the antiferromagnetic part of 
the superspin, is located at q = Q, where Q = (vr/a, tt/o) 
is the antiferromagnetic ordering vector. The uniform 
superconducting state (/ = 1) can be characterized by a 
superspin n = (/ cos (/), 0, 0, 0, /sin </>), and the tt mode is 
a fluctuation Sn{t) cx (0, 0, 0, /e'"''*, 0) around the static 
solution, where z has been chosen as an arbitrary direc- 
tion for Sm. In this case with / = 1 we have (5m (x e'""', 
i.e. a sharp peak at w = and q = Q. 

In the presence of an applied magnetic field, the super- 
conductor will be penetrated by flux quanta, each form- 
ing a vortex with a flux h/2e by which the complex su- 
perconducting order parameter tjj acquires a phase shift 
of 27r when moving around the vortex. In YBa2Cu307 
the vostices arrange themselves in a triangular vortex 
latticala with an area of the hexagonal unit cell given by 
A — h/2eB and consequently a lattice constant given by 
d — 37^'^y/h/eB. In the work by Arovas et Bruus 
et aL,ll3 and Alama et al&3 the problem of Abrikosay. 
vortices was studied within the SO (5) model of Zhang.EHl 
In the center of a vortex core, the superconducting part 
of the order parameter is forced to zero. This leaves two 
possibilities: i) either the vortex core is in a metallic 
normal state (as it is the case in conventional supercon- 
ductors) corresponding to a vanishing superspin or ii) 
the superspin remains intact but is rotated from the sua. 
perconducting sector into the antiferromagnetic sector .Ell 
The prediction of the possibility of antiferromagnetically 
ordered insulating vortex cores is thus quite novel and 
allows for a direct experimental test of the SO (5) the- 
ory. However, the antiferromagnetic ordering of vortices 
is according to our knowledge still to be confirmed exper- 
imentally. In this paper we report a different consequence 
of the S0(5) theory in neutron scattering experiments; 
we consider the tt mode in the presence of vortices and 
show that the peak at q = Q splits into a ring with a ra- 
dius n/d centered at q = Q where it has zero amplitude. 
Consequently the splitting scales with magnetic field B 
as ir/d oc Vb. 

We start by considering just one vortex, then gener- 
alize the result to a vortex lattice. To make our calcu- 
lations quantitative, we consider YBa2Cu307 for which 
a = 3.8 A, K ~ 84, and ^ ~ 16A for the lattice con- 
stant, the Ginzburg-Landau parameter, and the coher- 
ence length, respectively. The order parameter can be 
written in the formllfl 



without loss of generality. Static numerical solutions for 
/(r) and thereby also m(r) in the presence of a vortex 
are derived as described in Refs. Due to the high 

value of K the absolute value / of the superconducting 
order parameter ip increases from zero at the center of 
the vortex (r = 0) to its bulk value (/ = 1) at a distance 
of the order ^ from the center. The antiferromagnetic 
order parameter follows from / since m — yl — 

Eor the tt mode in the presence of a vortex, Bruus et 
aLt3 found that the fluctuation of the superspin is 



Sn{r,t) = (0,0,0,56'/(r)cos(/)re''^'*,0), 



(2) 



where the small angle 66 by which n rotates into the 
antiferromagnetic sector is undetermined. Since the ex- 
citation depends on / and not on m it is a de-localized 
excitation with zero amplitude at the center of the vor- 
tices and in terms of energy it actually corresponds to 
an energy at the bottom edge of the continuum of an 
effective potential associated to the vortices.tS 

For an isotropic spin space, the magnetic scattering 
cross-section for neutrons is proportional to the dynamic 
structure factor, which is the Fourier transform of the 
spin-spin correlation function (see e.g. Ref. pO| ), 

/oo 
dfe^'^*^e-''i-(^-^)(SR(i).SR.(0) 
RR' 

(3) 

To make a connection to the SO (5) calculations we 
make the semiclassical approximation (SR(i) • Sr'(O)) « 
(SR(t)).(SR.(0)) so that 
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-'i(q+Q)-(R-R') 



(4) 



where m(R, t) — e**^'^SR(i) is the antiferromagnetic or- 
der parameter which enters the superspin n. 

With a superspin given by n(r, t) = n(r) -I- 6n{r, t) 
the dynamical structure factor has two components — 
an elastic and an inelastic. The elastic component 



5ci(q,t^) 



R 



-l(q+Q)-R 



m(i?) 



2tt5{uj) 



(5) 



is located at q = Q and has a width ^ tt/^. In elastic 
neutron scattering experiments the observation of this 
peak would directly prove the antiferromagnetical order- 
ing in vortex cores. 

The inelastic contribution is 



n(r) = (/(r)cos(?!)r,0,m(r),0, /(r)sin(?!)r), (1) 

where 0r = arg(r). The isotropy of the antiferromagnetic 
subspace allows us to choose m to lie in the y-direction 



5in(q,c^) = {sey 



J2e-^^^+^^^fiR)cos^n 



R 



x2tt6{uj — lOtt). 



(6) 
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For q = Q the phase factor e~*('i+*^)^ vanishes, and the 
cosine factor makes the different terms in the summation 
cancel pairwise so that 5in(Q,W7r) = 0. The presence 
of a single vortex moves the intensity away from q = Q 
and a ring-shaped peak with radius Sq ^ tt/L centered at 
q = Q is formed, L ^ \J~A being the size of the sample. 
In the scmiclassical approximation the zero amplitude at 
q = Q is a topological feature, which is independent of 
the detailed radial form /(r) of the vortex. This robust- 
ness relies on the identification of the tt mode as being 
proportional to the superconducting order-parameter (in- 
cluding its phase). Quantum fluctuations may add some 
amplitude at q = Q, but such an analysis beyond leading 
order is outside the scope of this work. 

It is interesting to see how this result compares to 
predictions based on the BCS theory. The neutron 
scattering cross-section is given by the spin suscepti- 
bility, which for a homogeneous (vortex free) super- 
conductojrJias been calculated via the BCS-Lindhard 
functionJ3fl Here we briefly consider how the BCS co- 
herence factoij-ljUfcWfc+g — WfeMfe+g]^ appearing in the Lind- 
hard functions is modified by the presence of vortices. 
In a scmiclassical approximationEj the spatial variation 
of the superconducting phase 0(r) leads to a coher- 
ence factor of the form [Mfc(ri)e*"*(''i)/2wfc+^(r2)e''^('"^)/2 - 
?;fc(ri)e^'^(''i)/2^t^_^_^(i.2)e^'^('"2)/2]2. Therefore in contrast 
to Eq. (^ the superconducting phase does not separate 
in the two spatial positions, and consequently the spatial 
average in general is not zero at q = Q. It thus ap- 
pears that the above mentioned ring-shaped peak in the 
dynamic structure factor is special for the S0(5) model. 

We now generalize the single- vortex S0(5)-result to the 
case of a vortex lattice. For non-overlapping vortices we 
construct the full superconducting order parameter by 



V'(r) = /(r)e' 



i0(r) 



(7) 



where the Vj denote the positions of the vortices. The 
function /(r) — Hj /('^ ^ I'j) 1 except for close to 
the vortices where it dips to zero. Also the phase 
(^(r) = arg(r — r^) has by construction the period- 
icity of the vortex lattice (modulo 2-k) and the contour 
integral §^ dl-V0(r) equals 27rn where n is the number of 
vortices enclosed by the contour C. In the limit of non- 
overlapping vortices we can capture the main physics by 
considering the single vortex solution within a unit cell 
of the vortex lattice. We comment on the inclusion of 
the entire vortex lattice further on, but for now we re- 
strict the summation in Eq. (||) to lattice sites R inside 
the vortex lattice unit cell. In Fig. |l| we show the result 
for a magnetic field B = lOT. As seen, the presence of 
vortices moves the intensity away from q = Q and a ring 
shaped peak with radius 8q centered at q = Q is formed. 
We note that the only relevant length scale available is 
the vortex lattice constant d and consequently we expect 
that 8q = Ti/d. Since d = 'i'^^'^^Jh/eB we consequently 
expect that 8q = Z^I^-K^f^WJh ~ 0.008 x (tt/o) v/5/[T]. 




FIG. 1. Plot of the dynamic structure factor aX, lo = uj-n as 
a function of g along the (vr, 7r)-direction for B = 10 T. The 
inset shows the almost isotropic response in the g-plane with 
the arrow indicating the (tt, 7r)-direction. 



Had we included all the vortex lattice unit cells in our 
analysis, the structure factor of the hexagonal vortex lat- 
tice would have led to a breaking of the ring in Fig. Qinto 
six sub-peaks sitting on top of the ring. In a real experi- 
ment these sub-peaks could easily be smeared back into 
a ring-shaped scattering peak if either the vortex lattice 
were slightly imperfect or if the resolution of the spec- 
trometer were too low. To describe the main effect of the 
S0(5) theory we therefore continue to use the single unit 
cell approximation. 

In Fig. 1^ we show the splitting as a function of the 
magnetic field and indeed we find the expected scaling 
with a pre-factor confirming that the splitting is given 
by (5q = Tr/d. The full width half maximum of the ring is 
given by F ~ 3.1 x Sq = 3.1 x ir/d. 

In Fig. H we show the amplitude of the ring as a func- 
tion of magnetic field. The amplitude approximately de- 
creases as 1/B with the magnetic field, but with a small 
deviation. This deviation makes the q-integrated inten- 
sity, which is proportional to the amplitude times {Sq)^, 
decrease as I{B)/I{0) ~ 1 - 0.004 x B/[T] which re- 
fiects that the area occupied by vortices increases linearly 
with B and consequently the superconducting region de- 
creases linearly with B. In fact, the reduction is given by 
A~^27r Jrdrm^ir) ~ 0.004 x B/[T], where the integral 
gives the effective area of the vortex. The reduction in 
integrated intensity should be relatively easy to observe 
experimentally, but is not a unique feature of the S0(5) 
model. Thus, while it will aid to prove that the tt reso- 
nance only resides in the superconducting phase, it will 
not clearly distinguish between different theories. 

In order to discuss the experimental possibilities for 
testing our predictions, we note that the original obser- 
vation of the zero-field tt resonance was an experimen- 
tal achievement and hence that the experiment proposed 
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FIG. 2. Plot of the peak splitting 5q as a function of 
the magnetic field B. The calculated splitting (•) has the 
expected B^^^ behavior and the numerical pre-factor corre- 
sponds to a splitting 5q = n/d where d is the vortex distance. 



here constitutes a great challenge. However, since the 
first observation of the tt resonance in 1991, the field of 
neutron scattering has developed considerably. To ob- 
serve the ring-like shape (see inset of Fig. |^) of the exci- 
tation would require a resolution better than ir/d along 
two directions in reciprocal space, which seems unachiev- 
able with current spectrometers. However, the overall 
width of the ring can in fact be measured with good res- 
olution along just one direction in the reciprocal plane. 
Scans along this direction (as in Fig. ^ could then reveal 
a broadening of ^ 3.1 x ir/d. With a sufficiently opti- 
mized spectrometer we believe this to be possible, and 
the reward is a stringent test of a quantitative predie^ 
tion of the S0(5) theory. We note that Bourges et alx3 
have investigated the tt resonance in a magnetic field of 
B — 11.5 T and report a broadening in energy, but do 
not report data on the q-shape. 

In conclusion we have found that within the SO (5) 
model, the tt resonance splits into a ring centered at 
q = (tt/c, 7r/a) in the presence of a magnetic field. The 
ring has the radius n/d and full width half maximum of 
about 3.1 X Tr/d, where d is the vortex lattice constant. 
Consequently the splitting is found to scale with the mag- 
netic field as B^/^. We emphasize that the amplitude of 
the TT resonance is zero at q = (tt/o, n/a) in the presence 
of a magnetic field. 
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Ole R0mer Grant No. 9600548. 
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FIG. 3. The peak height plotted versus the magnetic field 
B. The calculations (•) almost fit a 1/B dependence. 
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